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I. INTRODUCTION 

The muon decay into an electron and a pair of neutrinos, /i — > ev^Pe, is a classical process in particle physics. 
Although the high energy frontier has moved up from the energy scales comparable to the muon mass, precision 
physics of muons remains an interesting and inspiring source of information about the Standard Model (SM) and its 
possible extensions 0. 

Among very different experiments with muons that include the measurement of muon anomalous magnetic moment, 
muon lifetime, fj, — ► e*y branching ratio, and the muon to electron conversion rate in muonic atoms, our attention 
will be focused on the TWIST experiment j^, |j| at TRIUMF, where the electron energy spectrum in muon decays 
is planned to be measured to determine Michel parameters with the precision of ~ 10~ 4 . In order to confront these 
measurements with the SM predictions and look for the signs of New Physics, one needs an adequately accurate 
calculation of the electron energy spectrum within the SM. 

Calculations of the electron energy spectrum in muon decay have a long and interesting history that dates back 
to the very early days of QED and the physics of weak interactions (see e.g. Ref. [Q for historical recollection). In 
spite of tremendous progress in precision calculations, the 0{a 2 ) radiative corrections to total muon lifetime have 
been computed only recently B , and the calculation of similar corrections to electron energy spectrum has even not 
been attempted. One reason for this is that, in contrast to total lifetime, the electron energy spectrum in muon decay 
can not be computed for the vanishing mass of the electron since terms enhanced by the large logarithm ln(m M /m e ) 
are present. These terms (excluding the ones that are related to the on-shell definition of the fine-structure constant 
commonly used in QED) cancel out in the calculation of the total rate which, for this reason, becomes somewhat 
simpler. 

At order 0{a 2 ), corrections to the spectrum have double-logarithmic O (In 2 (rri^/m 2 ,)) and single-logarithmic 
0(\n(m 2 ^/m 2 e )) enhanced terms and it is the purpose of this Letter to present the calculation of those. The 
double-logarithmic terms were computed recently in Ref. It was pointed out there that the single-logarithmic 
0(ln(m 2 /m 2 )) terms are required to match the precision of the TWIST experiment. Motivated by these considera- 
tions, we decided to perform this calculation. In order to accomplish this, we make use of perturbative fragmentation 
function approach borrowed to a large extent from QCD studies of heavy quark fragmentation in e + e~ collisions. 
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II. PRELIMINARIES 



According to QCD factorization theorems 0, which can be transformed to QED in a straightforward way, a 
differential cross section for producing a particle of a given type with a certain fraction of the initial energy can be 
written as a convolution of the hard scattering cross section computed with massless partons and the fragmentation 
function that describes the probability that a massless parton of a given type fragments to the physical particle in 
the final state that is being observed. 

If one considers the process in which a heavy quark (i.e. niQ ^> Aqcd) is produced with large energy and its 
energy is measured, one can identify the massive quark with the physical particle in the final state, in the sense of 
the preceding discussion. It has been shown in QCD that in this case the perturbative fragmentation function can be 
defined and that this function absorbs all the terms that are singular in the limit when the mass of the heavy quark 
goes to zero §, |, |(|. 

Specializing to the case of muon decay, and working along these ideas, it is possible to write down the formula for 
the electron energy spectrum in muon decay in the following way: 



dr ■c-^ f dz dTj fx , 

— {x,m f _ l ,m e )=}_^ / — — {z,m^^ f )VA -,fj. f ,m e ), (1) 



dx — / j 



where z = 2E/m fl is the fraction of energy carried away by a parton j in muon decay, x is the same quantity for the 
observed physical massive electron, dTj/dz is the energy distribution of the massless parton of type j in muon decay 
computed in the MS scheme, Dj is the fragmentation function of the parton j fragmenting into the massive electron, 
and fXf stands for the factorization scale. Note that the terms that are suppressed by the ratio of the electron to muon 
masses, m 2 /m 2 l , can not be described by Eq. (|l|). However, since these terms are known, both for Born and 0(a) 
corrected electron energy spectrum ]IT| , |l2| , Eq. (|l|) is quite adequate for anticipated level of experimental precision. 

As we mentioned earlier, partonic decay rate dTj/dz has to be computed in the MS scheme. This is important in 
that this requirement goes beyond standard ultraviolet renormalization, since dTj/dz is not finite because of collinear 
singularities. These singularities are removed from dTj/dz by conventional renormalization in the MS scheme and 
large collinear logarithms associated with them are absorbed into the fragmentation function Vj. 

The perturbative expansion for the energy distribution of the massless partons reads 

^ ^(*, m p> M/) = Af (z) + °MA? > K , f ,z) + (HM) * A? M/ , z), (2) 

where Tq — G|.m^/(967r 3 ), Aj(z) = z 2 (3 — 2z)8j e , a(/i/) is the MS renormalized fine structure constant, and we 
have neglected terms of order 0(a 3 ) and higher. The MS fine structure constant will be later converted into the 
on-shell coupling fine structure constant a ~ 1/137.036. 

Before giving explicit expressions for perturbative coefficients, we would like to describe a simple idea (previously 
exploited in QCD studies) that allows us to compute the a 2 \n.{m^/m e ) enhanced terms without doing explicit two- 
loop calculation. Since dTj/dz is computed for massless partons, it has just two energy scales, the mass of the muon 
m M and the factorization scale fj. f . This means that the only logarithms that arise there are the logarithms of the form 
h^m^/zz/). Therefore, by choosing [if ~ m M , we effectively eliminate large logarithms from perturbative coefficients 
of the energy spectrum computed with massless partons and move all the large logarithms to the fragmentation 
function T>j. On the other hand, the fact that the fragmentation function is process-independent and also satisfies the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation as far as its dependence on the factorization 
scale is concerned, allows us to use some formulas known from previous QCD studies and compute large logarithmic 
corrections to the electron energy spectrum at order 0{a 2 ) without performing explicit two-loop calculations. 

Let us consider the fragmentation function T>j(x, fif,m e ) which describes the probability that the massless parton 
j converts into the physical electron of the mass m e . This function satisfies the DGLAP evolution equation: 

l 

dVi(x,fif,m e ) m-^ f dz fx \ 

= J — P ii( Z > a W)) V A ->»f, m e h ( 3 ) 

where Pji stands for usual time-like splitting functions which, to the order we work to, can be written as 

P Tl {xMH)) - + (^) 'l$H*) + 0(aP) • (4) 
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Eq. (|^) can be solved as power series in a if the initial condition for the function T>j at some scale fig is provided. 
This initial condition can be obtained from QCD studies of the heavy quark fragmentation H which, when translated 
to QED, imply that the fragmentation of the massless electron into the physical electron is described by 



Zir 



ln-4-21n(l 

mi 



(5) 



and the fragmentation function of the photon fragmenting into the physical electron is 



£>' m (2,^0,TO e ) 



2tt 



(, 2 + (l-,) 2 )ln^. + 0( a 2 ) 



(6) 



As will be clear from the following discussion, the 0(a 2 ) terms in the initial conditions for fragmentation functions are 
not needed for our purposes since our choice of the initial scale /zo ~ m e guarantees that no large logarithms appear 
in the initial condition for T>j. One should also notice that at the lowest order O(a ) the fragmentation function does 
not contain large logarithm and, for this reason, the second order coefficient in dYj/dz is not needed as well. On the 
other hand, the order 0(a) coefficients in dTj/dz have to be known exactly and they read 



4 X) (*) 



A«(z) = (2z 2 (2z-3)ln 




1 

2z 2 (2z - 3) (4C 2 - 4Li 2 (z) + 2 In 2 z - 31nzln(l - z) - ln 2 (l - z)) 
5 „ _ , 34 



2z 



3 

5 23 

6 ~ T' 



13z 2 + —z 3 ) ln(l -z) + 



4z - 2z z - 6z 3 In z 
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—z 
2 



hi 



Tfl 



2z - 2z 2 + -z 3 
3 



Inz 



10 

t 



•4z 



60 

12* 



2z 2 



:z 3 . 



(7) 



(8) 



With these preliminary remarks out of the way, we can proceed with the computation of the fragmentation function 
and use it to calculate the electron energy spectrum in muon decay. 



III. FRAGMENTATION FUNCTION 



In this Section we discuss the computation of the fragmentation function. For this purpose, we need to solve the 
DGLAP equation (g) in a way consistent with the initial conditions. Since we solve this equation perturbatively, we 
can express the running fine structure constant in the MS scheme through the fine structure constant defined in the 
on-shell scheme, because the use of the on-shell renormalization scheme is the common practice when dealing with 
QED corrections. To order 0(a 2 ), the well known relation between the MS and the on-shell coupling constants reads 



a 2 (J-) 
a^ f ) = a + -ln—. 

Solving Eq. @ by iterations, we obtain 

V e (x,^ f ,m e ) = 5(1 - x) + ^- \LPf^{x) + d 1 (x,(M},m e fj 



(9) 



l -P£>®P£\x) 



V^x,n f ,m e ) = —LP^(x)+0(a 2 ) 



P^ ® di(ar) + P£\x)] ) + O {a 2 L°, a 3 ) 



(10) 
(11) 
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where L — h^/i 2 ://! 2 ,), and the convolution operation is defined in the standard way: 



A ® B(x) = J dz J dz' S(x - zz')A{z)B{z') = J y A(z)B (j- 

x 

We now give explicit expressions for the splitting functions used in Eqs. (lOfTl]). At leading order they read 



1 + x 2 
1-x 



l + (l-x) 2 



P^{x)=x 2 + {l-xf 



(12) 



(13) 



At next-to-leading order the time-like splitting functions have been derived for QCD in Refs. [H| [T(| 

and, by choosing appropriate color structures, can be translated to QED. Since, experimentally, one will probably 
distinguish events with one and more electrons in the final state, we decided to split the corresponding second order 
function P$ (x) into four parts in the same way as in Ref. |l7j : 



(14) 



Here Pel n \x) is determined by the set of Feynman diagrams with pure photonic corrections (i.e. no additional 
electrons in the final state or closed electron loops in virtual corrections); Pel' NS \x) is related to corrections due to 



non-singlet real and virtual e + e pairs; Pel' S \x) stands for the singlet pair production contribution; and Pel' lnt ^ 
describes the interference of the singlet and non-singlet pairs. These functions read 



(x) 



Pg«\x) = 8(1 -x) (--3& + 60 



1 + x 2 
1-x 



x) In 2 x -+ 


2x In x — 


3a; + 


2- 






1\ 20 


1 




2 1 + x 2 


6J ~ ~9~ 


1 - 


X 


, 3 l-i 



21nxln(l - x) - 2 In 2 a; - 2Li 2 (l - x) 



2 22 



Pg^{x) 

P£ S) (x) = (1 ------ .rihr.r -!- j - ').,- - j hl.r -S - ^- 4.r - — r 

Pg^Xx) 
where we have used 
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1 + x 2 
1-x 



7 



2Li 2 (l - i) + -lni - -(1 + x)lnx- 7 + 8x, 



OC ^ 



Li 2 (x) 



dz 



ln(l - z) 



(15) 
(16) 

(17) 
(18) 

(19) 



Finally, we give explicit formulas for various convolutions, which appear in Eq. (pi 



Pg><S>P£Hx)=6(l-x)[l-4&) + 



1 



1-x 



6 + 81n(l -x) 



1 -x 



lnx + (1 + x)[31nx - 41n(l - x)] - x - 5, 
P<° } ® P e ( ° } (x) = i— - (4 + 7x + 4x 2 ) + 2(1 + x) In x, 



3x 



(20) 
(21) 



p(°) ® dl (x) = pj°) ® pj°)(x) (in i£ - i) + «y(i - x) (| - 80 

7 + 8C2 -61n(l-a;) - 121n 2 (l - x) 



1-x 



1-x 



lnxln(l — x) 



+ (1 + x)[61n 2 (l - x) -61nxln(l - x) - 2Li 2 (l - x) - 4C 2 ] + 2xlnx 

+ (7-x)ln(l-x)-y-|z. (22) 
Using these results in Eqs.(nOUl^), one obtains explicit expressions for the fragmentation functions D en (x,[if,in e ). 
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IV. ELECTRON ENERGY SPECTRUM 

To obtain the electron energy spectrum we have to convolute the fragmentation functions Eqs. (JT^JTll) with dtj/dz, 
and we have given all the necessary ingredients to do that in the previous Sections. Before presenting corrections to 
the electron energy spectrum, let us note that the dependence on the factorization scale cancels out explicitly in the 
final result up to the terms that are not enhanced by any large logarithm and for this reason are beyond the scope of 
this paper. 

As we explained earlier, we split the final result into four different pieces: pure photonic radiative corrections, 
non-singlet pair radiative corrections, singlet pair radiative corrections, and corrections due to the interference of 
singlet and non-singlet pairs. Following this separation, we write the result for the electron energy spectrum in muon 
decay as 



J_^I = A(7) + A (NS) + A (S) + A (int)_ 

r da; 



(23) 



Let us start with pure photonic radiative corrections. Computing the convolutions, we arrive at the following result: 



2vr" 



2tt 



1 



(1,7)/ 



(24) 



where the dots stand for 0(a 2 ) terms without the logarithmic enhancement as well as higher order terms in the 
expansion of the electron energy spectrum in the fine structure constant. The O(a ) energy spectrum is given by 
fo(x) = x 2 (3 — 2x). The O(a) correction to it, fi(x), was calculated in Ref. [jll). The coefficient of the double-log 
term is given by 



f^\x) = ix 2 (3-2x) 



10 



iln 2 .x + ln 2 (l - x) - 21nxln(l - x) - Li 2 (l - x) - C2 



^- + %x - 16a: 2 + ^x 3 I ln(l - x) + I -~ - 2x + 8x" - ^x 3 I ma; 



3 

11 17 

1 ; 

36 6 



,.~ : *2 



G 



-a: 2 - —a: 3 



(25) 



and is therefore in agreement with recent results in Ref.||. The coefficient of the single-log enhanced term for pure 
photonic corrections, which is one of the new results in this Letter, reads 



2 (1 ' 7) (a;) = 2ir 2 (3 - 2x) ^-2Li 3 (a;) - 2Si, 2 (x) + 2Li 2 {x) ln(l - x) + 2U 2 {x) lnx 

+ 51na;ln 2 (l - x) - bin 2 x\n{\ - x) + 2 In 3 x - 2( 2 ln(l - x) - 2( 2 \nx + 7( 3 



Li 2 (x) 



10 



14x - AOx 2 



3 3 

+ ln 2 (l - x) (^-I2x - Ax 2 + 8x 3 J + In 2 x 

J 17 53 64 2 , - 
+ m(l - x)\ x H x 2 - 12x 3 



92 \ ( 25 92 
a; 3 ) +lnxln(l - x)[ h 32a: - 54a; 2 H : 



25 
12 



9 70 o 
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98 

+ C2 1 — 7T - 2a; + 35a; 2 - —a; 3 ) + — ^—x + ^-x z - ^x 



211 287 



83 n2 559 3 



216 12 



54 



44 



(26) 



where we used 



Li 3 (a;) 



dz 



U 2 (z) 



Si,2 (x) = - I dz 



ln 2 (l 



(27) 



The correction due to non-singlet electron-positron pairs, including effects of the running coupling constant, reads 



A ( NS ) 



ijf NS >(,)ln 2 ^<< 



/: 



(1,NS) 



(x) In 



(28) 
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with 



/ 2 (0 ' NS) (x) = 2x 2 (3-2x)ln 



1-x 



6 



2x - Ax 2 + -x 3 , 



/ 2 (1,NS) (x) = 2x 2 (3- 2x)[ -2Li 2 (l - x) - ~ lnxln(l - a;) + ^ln 2 (l - x) -ln 2 x- ~Ca 



,10 4 46 o „ , 

+ ln(l - x)[ x x 2 + 12a; 3 

1 9 3 3 



In a; I - + -x + 8x 2 
9 3 
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11 19 



100 2 64 , 
-x - — x 3 . 
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Next, we derive the result for the singlet pair correction. Writing 



(29) 



(30) 



A (S) 



(0 



1 "(°. S )^M 2f TO '\ , Al,S) t ^^( m l 



2/ --(x)ln^-f ) + /--(x)ln 



(31) 



we obtain 



/ 2 (i ' s) w 



2 17 „ 14 2 

1 h 3x x 

3x 9 3 



Li 2 (1 — x) + In x ln(l — x) 



9 V3 
5 
3 



4x + 4x 2 In x 



4x + 4x 2 ) + In 2 x I x + 6x + 4x i 



1 67 43 



77 



10 



3x 9 + 18 X+ 18* + 9 1 
Finally, for the interference term one finds 



, w'!7 2 o 14 2 8 3 \ , /8 4 5 19 2 



(32) 



(33) 



A (int) = / 

2tt 



/^'^(x)ln -£ 



(34) 



and 



/ 2 (1,int) (x) = 2x 2 (3-2x)^2Li 3 (l-x)-4Si i2 (l-x)-2Li 2 (l-x)lnx 
+ Li 2 (l-x)^+4x-26x 2 + yx 3 j + In 2 x f -9x 2 + yX 3 



lnx x 

3 3 



2<S 



62 41 



55 2 104 , 
x x H x . 



(35) 



V. CONCLUSIONS 



By applying techniques developed for perturbative QCD to Quantum Electrodynamics, we computed 0(a 2 ) correc- 
tion to the electron energy spectrum in unpolarized muon decay, keeping all the terms enhanced by the logarithms of 
the muon to electron mass ratio. The double logarithmic corrections 0(a 2 (In 2 m 2 /m 2 )) are in agreement with recent 
result in Ref. [||. The single logarithmic corrections 0(a 2 ln(m 2 /m 2 )), the new result of this Letter, are important 
to match the precision requirements of the TWIST experiment. To illustrate significance of single logarithmic terms, 
in Fig. |l| we plot both double and single 0(a 2 ) logarithmic corrections defined as 



4°' tot) (x) 



hr 



C (l,tot) 



(X) 



m (£) 2 [ /2(1,7)(x) + /2(1,NS)(x) + f2hS){x) + f2lMt){x) 



In 



(36) 
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as the function of the energy fraction carried away by electron in muon decay. 

As follows from Fig. [l], the O (a 2 \n{m 2 /m^)) corrections computed in this Letter are required for the theoretical 
prediction with the precision 1CP 4 , which is the benchmark precision for the TWIST experiment. Moreover, within 
the acceptance region of the TWIST experiment, 0.3 < x < 0.98, the magnitude of O {a 2 \n.(m fl /m e )) corrections is 
in fact comparable to the magnitude of (naively the largest) O (a 2 ln 2 (m A1 /m e )) terms. 
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FIG. 1: Double and single logarithmic corrections as a function of x. 

The fact that sub-leading logarithmically enhanced corrections are larger than the precision of the TWIST experi- 
ment, seems to indicate that full calculation of the 0(a 2 ) corrections to the electron energy spectrum is very desirable 
and that without such a calculation an intrinsic theory uncertainty can not be pushed below few x 10~ 4 . 

Complete calculation of 0(a 2 ) correction to the electron energy spectrum in muon decay is a very difficult task and 
it remains to be seen if it can be done. In spite of that, it is possible to extend the analysis of this Letter in a number 
of ways to further improve on the theory prediction. First of all, it is straightforward to compute the electron energy 
spectrum in polarized muon decay, by using techniques described in this Letter. Furthermore, O (a 3 In 3 {rn^/m 2 )) 
corrections can be obtained from the DGLAP equation. Also, resummation of the corrections that are singular in the 
limit x — > 1 can be performed and the influence of this resummation on theoretical predictions for the spectrum can 
be explored. These studies as well as detailed discussion of the present theoretical uncertainty in the electron energy 
spectrum in polarized muon decay will be presented elsewhere. 
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